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Abstract 



We present an invariant definition of the hypoelliptic Laplacian on sub-Riemannian 
structures with constant growth vector using the Popp's volume form introduced by 
Montgomery. This definition generalizes the one of the Laplace-Beltrami operator 
in Riemannian geometry. In the case of left-invariant problems on unimodular Lie 
groups we prove that it coincides with the usual sum of squares. 

We then extend a method (first used by Hulanicki on the Heisenberg group) to 
compute explicitly the kernel of the hypoelliptic heat equation on any unimodular 
Lie group of type I. The main tool is the noncommutative Fourier transform. We 
then study some relevant cases: SU(2), SO(3), SL(2) (with the metrics inherited 
by the Killing form), and the group SE(2) of rototranslations of the plane. 
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1 Introduction 

The study of the properties of the heat kernel in a sub-Riemannian manifold drew an increasing 
attention, since the pioneer work of Hormander [28]. 

Since then, many estimates and properties of the kernel in terms of the sub-Riemannian dis- 
tance have been provided (see [9,11,7,13,20,31,35,42,47] and references therein). 
In most of the cases the hypoelliptic Laplacian appearing in the heat equation is the the sum 
of squares of the vector field forming an orthonormal frame for the sub-Riemannian structure. 
In other cases it is built as the divergence of the horizontal gradient, where the divergence is 
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defined using any C°° volume form on the manifold (see for instance [46]). 
The Laplacians obtained in these ways are not intrinsic in the sense that they do not depend 
only on the sub-Riemannian distance. Indeed, when the Laplacian is built as the sum of squares, 
it depends on the choice of the orthonormal frame, while, when it is defined as divergence of 
the horizontal gradient, it depends on the choice of the volume form. 

The first question we address in this paper is the definition of an invariant hypoelliptic Lapla- 
cian. To our knowledge, the first time in which this question was pointed out was in a paper 
by Brockett [8]. Many details can be found in Montgomery's book [38]. 

To define the intrinsic hypoelliptic Laplacian, we proceed as in Riemannian geometry. In Rie- 
mannian geometry the invariant Laplacian (called the Laplace-Beltrami operator) is defined as 
the divergence of the gradient where the gradient is obtained via the Riemannian metric and 
the divergence via the Riemannian volume form. 

In sub-Riemannian geometry, we define the invariant hypoelliptic Laplacian as the divergence 
of the horizontal gradient. The horizontal gradient of a function is the natural generalization 
of the gradient in Riemannian geometry and it is a vector field belonging to the distribution. 
The divergence is computed with respect to the sub-Riemannian volume form, that can be de- 
fined for every sub-Riemannian structure with constant growth vector. This definition depends 
only on the sub-Riemannian structure. The sub-Riemannian volume form, called the Popp's 
measure, was first introduced in Montgomery's book [38], where its relation with the Hausdorff 
measure is also discussed. The definition of the sub-Riemannian volume form is simple in the 
3D contact case, and a bit more delicate in general. 

We then prove that for the wide class of unimodular Lie groups (i.e. the groups where the 
right- and left-Haar measures coincide) the hypoelliptic Laplacian is the sum of squares for 
any choice of a left-invariant orthonormal base. We recall that all compact and all nilpotent 
Lie groups are unimodular. 
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In the second part of the paper, we present a method to compute explicitly the kernel of the 
hypoelliptic heat equation on a wide class of left-invariant sub-Riemannian structures on Lie 
groups. We then apply this method to the most important 3D Lie groups: SU(2), SO (3), and 
SL{2) with the metric defined by the Killing form, the Heisenberg group H2, and the group 
of rototranslations of the plane SE(2). These groups are unimodular, hence the hypoelliptic 
Laplacian is the sum of squares. The interest in studying SU(2), 50(3) and SL(2) comes from 
some recent results of the authors. Indeed, in [10] the complete description of the cut and 
conjugate loci for these groups was obtained. These results, together with those of this paper, 
open new perspectives in the direction of clarifying the relation between the presence of the 
cut locus and the properties of the heat kernel in the spirit of the result of Neel and Strook 
[39] in Riemannian geometry. Up to now the only case in which both the cut locus and the 
heat kernel were known explicitly was the Heisenberg group [21,22,29]. 1 

The interest in the hypoelliptic heat kernel on SE(2) comes from a model of human vision. It 
was recognized in [15,41] that the visual cortex VI solves a nonisotropic diffusion problem on 
the group SE(2) while reconstructing a partially hidden or corrupted image. The study of the 
cut locus on SE(2) is a work in progress. Preliminary results can be found in [37]. 
The method is based upon the generalized (noncommutative) Fourier transform (GFT, for 
short), that permits to disintegrate 2 a function from a Lie group G to 1 on its components on 
(the class of) non-equivalent unitary irreducible representations of G. This technique permits 



1 The Heisenberg group is in a sense a very degenerate example. For instance, in this case the cut locus 
coincides globally with the conjugate locus (set of points where geodesies lose local optimality) and 
many properties that one expects to be distinct for more generic situations cannot be distinguished. 
The application of our method to the Heisenberg group H2 provides in a few lines the Gaveau- 
Hulanicki formula [21,29]. 

2 One could also say decompose (possibly continuously). 
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to transform the hypoelliptic heat equation into an equation in the dual of the group 3 , that 
is particularly simple since the GFT disintegrate the right-regular representations and the 
hypoelliptic Laplacian is built with left-invariant vector fields (to which a one parameter group 
of right-translations is associated). 

Unless we are in the abelian case, the dual of a Lie group in general is not a group. In the 
compact case it is a so called Tannaka category [26,27] and it is a discrete set. In the nilpotent 
case it has the structure of W 1 for some n. In the general case it can have a quite complicated 
structure. However, under certain hypotheses (see Section 3), it is a measure space if endowed 
with the so called Plancherel measure. Roughly speaking, the GFT is an isometry between 
L 2 (G, C) (the set of complex- valued square integrable functions over G, with respect to the 
Haar measure) and the set of Hilbert-Schmidt operators with respect to the Plancherel measure. 
The difficulties of applying our method in specific cases rely mostly on two points: 

i) computing the tools for the GFT, i.e. the non-equivalent irreducible representations of 
the group and the Plancherel measure. This is a difficult problem in general: however, for 
certain classes of Lie groups there are suitable techniques (for instance the Kirillov orbit 
method for nilpotent Lie groups [33], or methods for semidirect products). For the groups 
disccussed in this paper, the sets of non-equivalent irreducible representations (and hence 
the GFT) are well known (see for instance [43]); 

ii) finding the spectrum of an operator (the GFT of the hypoelliptic Laplacian). Depending 
on the structure of the group and on its dimension, this problem gives rise to a matrix 
equation, an ODE or a PDE. 

Then one can express the kernel of the hypoelliptic heat equation in terms of eigenfunctions of 
the GFT of the hypoelliptic Laplacian, or in terms of the kernel of the transformed equation. 
For the cases treated in this paper, we have the following (the symbol II means disjoint union): 

3 In this paper, by the dual of the group, we mean the support of the Plancherel measure on the set 
of non-equivalent unitary irreducible representations of G; we thus ignore the singular representations. 
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Group 


Dual of 
the group 


GFT 

of the hypoelliptic Laplacian 


Eigenfunctions of the GFT 
of the hypoelliptic Laplacian 


H 2 


R 


d 2 

— — - — X 2 x 2 (quantum Harmonic oscillator) 

dx l 


Hermite polynomials 


SU{2) 


N 


Linear finite dimensional operator related to the 
quantum angular momentum 


Complex homogeneous 
polynomials in two variables 


SO{3) 


N 


Linear finite dimensional operator related to 
orbital quantum angular momentum 


Spherical harmonics 


SL{2) 


R + II R + 
UN II N 


Continuous: Linear operator on analytic functions 

with domain {\x\ — 1} C C 
Discrete: Linear operator on analytic functions 

with domain { \x\ < 1} C C 


Complex monomials 


SE{2) 


R+ 


d 2 

— — — A 2 cos 2 (#) (Mathieu's equation) 
do 2 


27r-periodic Mathieu functions 



The idea of using the GFT to compute the hypoelliptic heat kernel is not new: it was already 
used on the Heisenberg group in [29] at the same time as the Gaveau formula was published 
in [21], and on all step 2 nilpotent Lie groups in [16,5]. See also the related work [34]. 



The structure of the paper is the following: in Section 2 we recall some basic definitions from 
sub-Riemannian geometry and we construct the sub-Riemannian volume form. We then give 
the definition of the hypoelliptic Laplacian on a regular sub-Riemannian manifold, and we 
show that the hypothesis of regularity cannot be dropped in general. For this purpose, we show 
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that the invariant hypoelliptic Laplacian defined on the Martinet sub-Riemannian structure 
is singular. We then move to left-invariant sub-Riemannian structures on Lie groups and we 
show that a Lie group is unimodular if and only if the invariant hypoelliptic Laplacian is the 
sum of squares. We also provide an example of a 3D non-unimodular Lie group for which the 
invariant hypoelliptic Laplacian is not the sum of squares. The section ends with the proof 
that the invariant hypoelliptic Laplacian can be expressed as 

m 

A sr = -^2L* x .L Xi , 

2 = 1 

where the formal adjoint L* x . is built with the sub-Riemannian volume form, providing a 
connection with existing literature (see e.g. [30]). The invariant hypoelliptic Laplacian is then 
the sum of squares when L x . are skew-adjoint 4 . 

In Section 3 we recall basic tools of the GFT and we describe our general method to compute 
the heat kernel of the hypoelliptic Laplacian on unimodular Lie groups of type I. We provide 
two useful formulas, one in the case where the GFT of the hypoelliptic Laplacian has discrete 
spectrum, and the other in the case where the GFT of the hypoelliptic heat equation admits 
a kernel. 

In Section 4 we apply our method to compute the kernel on H 2 , SU(2), SO (3), SL(2) and 
SE(2). For the Heisenberg group we use the formula involving the kernel of the transformed 
equation (the Mehler kernel). For the other groups we use the formula in terms of eigenvalues 
and eigenvectors of the GFT of the hypoelliptic Laplacian. 

The application of our method to higher dimensional sub-Riemannian problems and in par- 

4 This point of view permits to give an alternative proof of the fact that the invariant hypoelliptic 
Laplacian on left-invariant structures on unimodular Lie groups is the sum of squares. As a matter of 
fact, left-invariant vector fields are formally skew-adjoint with respect to the right-Haar measure. On 
Lie groups the invariant volume form is left-invariant, hence is proportional to the left-Haar measure, 
and is in turn proportional to the right-Haar measure on unimodular groups. 
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ticular to the nilpotent Lie groups (2, 3, 4) (the Engel group) and (2, 3, 5) is the subject of a 
forthcoming paper. 

2 The hypoelliptic Laplacian 

In this Section we give a definition of the hypoelliptic Laplacian A sr on a regular sub-Riemannian 
manifold M. 

2.1 Sub-Riemannian manifolds 

We start recalling the definition of sub-Riemannian manifold. 

Definition 1 A (n,m) -sub-Riemannian manifold is a triple (M, A,g), where 

• M is a connected smooth manifold of dimension n; 

• A is a smooth distribution of constant rank m < n satisfying the Hormander condition, 
i.e. A is a smooth map that associates to q G M a m-dim subspace A(q) ofT q M and V q G M 
we have 

spanffXi, [. . . [X k _ u X k ] . . ]}(q) | X t G Vec H (M)} = T q M (1) 
where Vec#(M) denotes the set of horizontal smooth vector fields on M, i.e. 

VeCtf(M) = {X G Vec(M) | X(p) G A(p) VpG M} . 

• g q is a Riemannian metric on A.(q), that is smooth as function of q. 

When M is an orientable manifold, we say that the sub-Riemannian manifold is orientable. 
Remark 2 Usually sub-Riemannian manifolds are defined with m < n. In our definition we 
decided to include the Riemannian case m = n, since all our results hold in that case. Notice 
that if m = n then condition (1) is automatically satisfied. 

A Lipschitz continuous curve 7 : [0,T] — > M is said to be horizontal if j(t) G A(7(t)) for 
almost every t G [0, T]. 
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Given an horizontal curve 7 : [0, T] — > M, the length 0/7 is 



*(7) = / o <** 



(2) 



The distance induced by the sub-Riemannian structure on M is the function 



d(qo,qi) 



inf{/(7) I 7(0) = qoy'jiT) = 91,7 horizontal}. 



(3) 



The hypothesis of connectedness of M and the Hormander condition guarantee the finiteness 
and the continuity of d(-, •) with respect to the topology of M (Chow's Theorem, see for 
instance [2]). The function d(-, •) is called the Carnot-Charateodory distance and gives to M 
the structure of metric space (see [6,23]). 

It is a standard fact that £(7) is invariant under reparameterization of the curve 7. Moreover, 
if an admissible curve 7 minimizes the so-called energy functional 



with T fixed (and fixed initial and final point), then v = y&y(t){j{t),j{t)) is constant and 7 
is also a minimizer of /(■). On the other side, a minimizer 7 of /(■) such that v is constant is a 
minimizer of E(-) with T = l{^f)/v. 

A geodesic for the sub-Riemannian manifold is a curve 7 : [0, T] — > M such that for every 
sufficiently small interval t 2 ] C [0, T], 7| [t t j is a minimizer of E(-). A geodesic for which 
g 7 (t)(7(t), j{t)) is (constantly) equal to one is said to be parameterized by arclength. 
Locally, the pair (A, g) can be given by assigning a set of m smooth vector fields spanning A 
and that are orthonormal for g, i.e. 



In this case, the set {X±, . . . , X m } is called a local orthonormal frame for the sub-Riemannian 
structure. When (A, g) can be defined as in (4) by m vector fields defined globally, we say that 
the sub-Riemannian manifold is trivializable. 





k(q) = span {Xi(q), ...,X m (q)}, g q (Xi(q), Xj(q)) = c%. 



(4) 
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Given a (n, m)- trivializable sub-Riemannian manifold, the problem of finding a curve mini- 
mizing the energy between two fixed points qo, q\ G M is naturally formulated as the optimal 
control problem 



It is a standard fact that this optimal control problem is equivalent to the minimum time 

problem with controls Ui,...,u m satisfying U\{t) 2 + . . . + u m (t) 2 < 1 in [0, T}. 

When the manifold is analytic and the orthonormal frame can be assigned through m analytic 

vector fields, we say that the sub-Riemannian manifold is analytic. 

We end this section with the definition of the small flag of the distribution ▲: 

Definition 3 Let A be a distribution and define through the recursive formula 



where A n+1 (q ) := A n (g ) + [A„(g ), A(g )] = 

= {X 1 (q ) + [X 2 ,X 3 ](go) | Xx(q),X 2 (q) G A n (q), X 3 (q) G A(q) V q G M} . The small flag of 
▲ is the sequence 



A sub-Riemannian manifold is said to be regular if for each n = 1,2,... the dimension of 
A n (<jo) = {/(<Zo) I f(q) £ A n (g) V q G M} does not depend on the point go £ M. 
A 3D sub-Riemannian manifold is said to be a 3D contact manifold if ▲ has dimension 2 
and A 2 (go) — T qo M for any point q G M. 

In this paper we always deal with regular sub-Riemannian manifolds. 
2.1.1 Left-invariant sub-Riemannian manifolds 

In this section we present a natural sub-Riemannian structure that can be defined on Lie 
groups. All along the paper, we use the notation for Lie groups of matrices. For general Lie 





9(0) = 9o, q(T) = q 1 . 



Ai:=A A n+1 := A„ + [A n , A] 



Ai C A 2 C . . . C A n C . . . 
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groups, by gv with g G G and v G L, we mean (L g )*(v) where L 5 is the left-translation of the 
group. 

Definition 4 Let G be a Lie group with Lie algebra L and p C L a subspace of L satisfying 
the Lie bracket generating condition 



Endow p wift a positive definite quadratic form (.,.). Define a sub-Riemannian structure on 
G as follows: 

• the distribution is the left-invariant distribution A(g) '■— gp; 



In this case we say that (G, A,g) is a left-invariant sub-Riemannian manifold. 

Remark 5 Observe that all left-invariant manifolds (G, A,g) are regular. 

In the following we define a left-invariant sub-Riemannian manifold choosing a set of m vectors 

{pi, . . . ,p m } being an orthonormal basis for the subspace p C L with respect to the metric 

defined in Definition 4, i.e. p = span {pi, . . . ,p m } and (pi,Pj) = We thus have A(g) = 

gp = span {gpi, . . . , gp m } and g g (gpi,gpj) = o~ij- Hence every left- invariant sub-Riemannian 

manifold is trivializable. 

The problem of finding the minimal energy between the identity and a point g\ G G in fixed 
time T becomes the left-invariant optimal control problem 



Lie p := span{[pi, [p 2 , \p n -i,Pn]]] \ Pi G p} = L. 



• the quadratic form g on ▲ is given by g g (w 1 ,f 2 ) := (g x Vi,g l v 2 ). 







g(0) = Id, g(T) = gi . 



Remark 6 This problem admits a solution, see for instance Chapter 5 of [12]. 
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2.2 Definition of the hypoelliptic Laplacian on a sub-Riemannian manifold 



In this section we define the intrinsic hypoelliptic Laplacian on a regular orientable sub- 
Riemannian manifold (M, A, g). This definition generalizes the one of the Laplace-Beltrami op- 
erator on an orientable Riemannian manifold, that is A0 := div grad 0, where grad is the unique 
operator from C°°(M) to Vec (M) satisfying g 9 (grad (f)(q),v) = d(p q {v) V q G M, v G T q M, 
and the divergence of a vector field X is the unique function satisfying divX/z = LxH where [i 
is the Riemannian volume form. 

We first define the sub-Riemannian gradient of a function, that is an horizontal vector field. 
Definition 7 Let (M, A,g) be a sub-Riemannian manifold: the horizontal gradient is the 

unique operator grad sr from C°°(M) to Vec#(M) satisfying g g (grad sr 0(g), v) = d<p q {v) V q G 
M, v G A(q). 

One can easily check that if {Xi, . . . X m } is a local orthonormal frame for (M, A,g), then 



The question of defining a sub-Riemannian volume form is more delicate. We start by consid- 
ering the 3D contact case. 

Proposition 8 Let (M, A, g) be an orientable 3D contact sub-Riemannian structure and {Xi, X 2 } 
a local orthonormal frame. Let X 3 = [Xl,^] anddX\, dX 2 , dX 3 the dual basis, i.e. dXi(Xj) = 
dij. Then fi sr := dX\ A dX 2 A dX 3 is an intrinsic volume form, i.e. it is invariant for a orien- 
tation preserving change of orthonormal frame. 



Proof. Consider two different orthonormal frames with the same orientation {Xi,X 2 } and 
{Fi,F 2 }- We have to prove that dX 1 A dX 2 A dX 3 = dY x A dY 2 A dY 3 with X 3 = [X 1 ,X 2 ], 
Y 3 = [Y U Y 2 \. We have 



grad sr = IXi (LxiflXi. 



( \ 

Y x 



( 



cos(/(g)) sin(/(g)) 



X 1 




V 



sin(/(g)) cos(/(g)) 
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for some real- valued smooth function /. A direct computation shows that 
Y 3 — X 3 + fiXi + J2X2 



(7) 



where fi and /2 are two smooth functions depending on /. 

We first prove that dX\ A dX 2 = dY\ A dYi. Since the change of variables {Xi, X2} 1— > {Y±, Y2} 
is norm-preserving, we have dXi A dX 2 (v, w) = dY\ A dY 2 (v,w) when v, w G A. Consider now 
any vector v = V\Xi + v 2 X 2 + = v[Yi + u 2 Y2 + & s a consequence of (7), we have 
V3 = v' 3 . Take another vector w = wiXi + w 2 X 2 + = w^Yi + io 2 ^2 + ^3^3 and compute 

dXiAdX 2 (v,w) = dXxhdX^v-VzXz^w-WzXz) = dYiAdY2(v—v s X 3 ,w—w 3 X 3 ) = dYxAdY^v.w), 

because the vectors v — v 3 X 3 , w — W3X3 are horizontal. Hence the two 2-forms coincide. 
From (7) we also have dY 3 = dX 3 + f[dX\ + f 2 dX 2 for some smooth functions f[, f 2 . Hence we 
have dYiAdY 2 AdY 3 = dX x AdX 2 AdY 3 = dX 1 AdX 2 A{dX 3 + f[dX 1 + f 2 [dX 2 ) = dX Y AdX 2 AdX 3 , 
where the last identity is a consequence of skew-symmetry of differential forms. □ 
Remark 9 Indeed, even if in the 3D contact case there is no scalar product in T q M, it is 
possible to define a natural volume form, since on ▲ the scalar product is defined by g and 
formula (7) guarantees the existence of a natural scalar product in (A + [A, A])/A. 
The previous result generalizes to any regular orientable sub-Riemannian structure, as pre- 
sented below. 

2.2.1 Definition of the intrinsic volume form 

Let = Eq C Ei C . . . C Ek = E be a filtration of an n-dimensional vector space E. Let 
ei, . . . , e n be a basis of E such that Ei = span {ei, . . . , e ni }. Obviously, the wedge product 
t\ A . . . A e n depends only on the residue classes 
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where < j < ra^+i, j = 1, . . . ,n. This property induces a natural (ie. independent on the 
choice of the basis) isomorphism of 1-dimensional spaces: 

n / k \ 

A n ^ = A ©(*W-i) . 



\t=l 



Now consider the filtration 



C Ai(g) C . . . C A fc (g) = T 9 M, dim A*(g) = n*. 

Let Xi, . . . , Xi be smooth sections of A = Ai; then the vector 

([X u [. . . ,Xi] . . .}(q) + A t _i(g)) G A^/A^g) 

depends only on -Xi(g) <g> . . . <g> Xi(q). 
We thus obtain a well-defined surjective linear map 

A{q) w - A,(g)/A^(g) 

A : 

Xi (g)®... i ([X 1 ,[...,X]...](g) + A 4 _ 1 (g); 

The Euclidean structure on A(g) induces an Euclidean structure on A(g)®* by the standard 

formula: 

(£1 ® ■ ■ ■ ® &> »7i ® • • • ® »?») = (fi, »7i) • ■ ■ (&, »/•), €j, Vj e A(g), j = 1, . . . , i. 
Then the formula: 

M = min{|el : £ G « e A^/A^g) 

defines an Euclidean norm on A i (g)/A i ~ 1 (g). 

Let be the volume form on A*(g)/A (g) associated with the Euclidean structure: 



i 

(u u vi A . . . v mt ) = det m i {(vj, vj')}jj/=i 
where m ; = n,j — = dim(A*(g)/A l_1 (g)). 
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Finally, the intrinsic volume form /i sr on T q M is the image of V\ A . . . A i/fe under the natural 
isomorphism 



Remark 10 The construction given above appeared for the first time in the book of Mongomery 
[38, Section 10.5]. He called the measure /i sr the Popp's measure. He also observed that a sub- 
Riemannian volume form was the only missing ingredient to get an intrinsic definition of 
hypoelliptic Laplacian. 5 

Once the volume form is defined, the divergence of a vector field X is defined as in Riemannian 
geometry, i.e. it is the function d\v sr X satisfying div sr Xfi sr = L x ^ S r- We are now ready to 
define the intrinsic hypoelliptic Laplacian. 

Definition 11 Let (M, A,g) be an orientable regular sub-Riemannian manifold. Then the in- 
trinsic hypoelliptic Laplacian is A sr <p := div sr grad sr 0. 

Consider now an orientable regular sub-Riemannian structure (M, A,g) and let {Xi, . . .X m } 
be a local orthonormal frame. We want to find an explicit expression for the operator A sr . 
If n = m then A sr is the Laplace Beltrami operator. Otherwise consider n — m vector fields 
X m+ i, . . . , X n such that {X^q), . . . , X m (q), X m+1 (g), . . . , X n (q)} is a basis of T q M for all q in 
a certain open set U . The volume form /i sr is fi sr = f\q)dXi A ... A dX n , with dXi dual basis 
of X\, . . . , X n : then we can find other n — m vector fields, that we still call X m+ i, . . . , X n , for 
which we have /i sr = dXi A ... A dX n . 

Recall that A sr satisfies (A sr <p)ii sr = Lxf^sr with X = grad sr 0. We have 



5 Montgomery did not use Popp's measure to get the intrinsic definition of the hypoelliptic Laplacian 
since there are, a priori, two natural measures on a regular sub-Riemannian manifold: the Popp's mea- 
sure, and the Hausdorff measure (see [36,40]). However, for left-invariant sub-Riemannian manifolds, 
both of them are proportional to the left-Haar measure. See Remark 16. 
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7(1 



Lx^sr = [d ((#, Xi)) A dXi A ... A dXi A ... A 



i=l 



+ X 4 )d (dXi A ... A dXi A ... A dX n ) } . 



Applying standard results of differential calculus, we have d ({d<f), Xi)) A dX\ A ... A dXi A ... A 



dX n = {-l) i+1 L Xi <j) /i sr and d (dX l A ... A dXi A ... A dX n ) = (-l) i+1 Tr (adX;) ^ sr , where 



the adjoint map is 



Vec ([/) 



Vec (U) 



adXi : 



X ^ [X 4 ,X]. 

and by Tr (ad Xj) we mean Z)?=i dXj([Xj, Xj]). Finally, we find the expression 



Til 



A sr <f) = £ (px^ + L x ^ Tr (adX t )) . 



(8) 



i=i 



Notice that the formula depends on the choice of the vector fields X m +i, . . . , X n . 

The hypoellipticity of A sr (i.e. given U C M and : U — * C such that A sr £ C°°, then is 

follows from the Hormander Theorem (see [28]): 
Theorem 12 Let L be a differential operator on a manifold M, that locally in a neighbor- 
hood U is written as L = Y^iLiL Xi + Lx , where X ,Xi . . . ,X m are C°° vector fields. If 
Lie g {X , Xi, . . . , X m } = T q M for all q EU , then L is hypoelliptic. 

Indeed, A sr is written locally as Y^Li L 2 Xi +Lx with the first-order term Lx = Y^Li Tr (adX) Lx v 
Moreover by hypothesis we have that Lie g {Xi, . . . , X m } = T q M, hence the Hormander theorem 
applies. 

Remark 13 Notice that in the Riemannian case, i.e. for m = n, A sr coincides with the Laplace- 
Beltrami operator. 

Remark 14 The hypothesis that the sub-Riemannian manifold is regular is crucial for the 
construction of the invariant volume form. For instance for the Martinet metric on IR 3 , that 

2 

is the sub-Riemannian structure for which L\ = d x + \d z and L 2 = d y form an orthonormal 
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base, one gets on M 3 \ {y = 0} 



A sr = (LO 2 + (L 2 ) 2 - -L 2 . 



y 



This is not surprising at all. As a matter of fact, even the Laplace-Beltrami operator is singular 
in almost- Riemannian geometry (see [3] and references therein). For instance, for the Grushin 
metric on R 2 , that is the singular Riemannian structure for which L\ = d x and L 2 = xd y form 
an orthonormal frame, one gets on IR 2 \ {x = 0} 



2.3 The hypoelliptic Laplacian on Lie groups 

In the case of left-invariant sub-Riemannian manifolds, there is an intrinsic global expression 
of A sr . 

Corollary 15 Let (G, A,g) be a left-invariant sub-Riemannian manifold generated by the or- 
thonormal basis {pi, . . . ,p m } C L. Then the hypoelliptic Laplacian is 



i=i 

where Lx i is the Lie derivative w.r.t. the field Xi = gpi. 

Proof. If m < n, we can find n — m vectors {p m +i, ■ ■ ■ ,p n } such that {pi, . . . ,p n } is a basis 
for L. Choose the fields Xi := gpi and follow the computation given above: we find formula 
(9). In this case the adjoint map is intrinsically defined and the trace does not depend on the 
choice of X m+1 , . . . , X n . □ 

The formula above reduces to the sum of squares in the wide class of unimodular Lie groups. 
We recall that on a Lie group of dimension n, there always exist a left-invariant n-form /i^ 
and a right-invariant n-form u# (called respectively left- and right-Haar measures), that are 




in 




(9) 
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unique up to a multiplicative constant. These forms have the properties that 
/ f{ag)n L (g)= / f(g)fi L (g), / f(ga)fi R {g) = / f(g)n R (g), for every / e L\G, R) and a e 

iG JG JG 

where L 1 is intended with respect to the left-Haar measure in the first identity and with respect 
to the right-Haar measure in the second one. The group is called unimodular if \xl and \i R are 
proportional. 

Remark 16 Notice that for left-invariant sub-Riemannian manifolds the intrinsic volume form 
and the Hausdorff measure fijj are left-invariant, hence they are proportional to the the left 
Haar measure [Ml- On unimodular Lie groups one can assume fi sr = [Ml = Hr = a^H, where 
a > is a constant that is unknown even for the simplest among the genuine sub-Riemannian 
structures, i.e. the Heisenberg group. 

Proposition 17 Under the hypotheses of Corollary 15, if G is unimodular then A sr (f) = 

Proof. Consider the modular function that is the unique function such that f G f{h~ 1 g)fin(g) = 
^{h) Jg f (9)1^11(9) f° r a ^ / measurable. It is well known that ^f(g) = det(Ad 5 ) and that 
ty(g) = 1 if and only if G is unimodular. 

Consider a curve "f(t) such that 7 exists for t — t : then j(t) = g e('~'°)' ?+o( * _ * ) with g — 7 (to) 
and for some 77 G L. We have 



^ det(Ad 7W )=Tr(^(Ad 90 ; 



d 



9Qe s>7 + o(s) 



(10) 



Tr (^Ad 9 -iAd go ad^j = Tr (ad,,) 



Now choose the curve j(t) = g§e tVi and observe that det(Ad 7 ( t )) = 1, then Tr(ad Pi ) = 0. The 
conclusion follows from (9). □ 

All the groups treated in this paper, (i.e. H 2 , SU(2), SO (3), SL{2) and SE{2)) are uni- 
modular. Hence the invariant hypoelliptic Laplacian is the sum of squares. A kind of inverse 
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result holds: 

Proposition 18 Let (G, A,g) be a left-invariant sub-Riemannian manifold generated by the 
orthonormal basis {pi, . . . ,p m } C L. If the hypoelliptic Laplacian satisfies A sr (p = Y%Li L\.<\), 
then G is unimodular. 

Proof. We start observing that A sr <j) = Y^iL\L 2 x .<j) if and only if Tr(ad Pi ) = for all % = 
1, . . . ,m. 

Fix g £ G: due to Lie bracket generating condition, the control system (6) is controllable, then 
there exists a choice of piecewise constant controls Ui : [0, T] — > R such that the corresponding 
solution 7(.) is an horizontal curve steering Id to g. Then 7 is defined for all t £ [0,T] except 
for a finite set E of switching times. 

Consider now the modular function along 7, i.e. ^(7^)), that is a continuous function, differ- 
entiable for all t £ [0, T]\E. We compute its derivative using (10): we have 4, det(Ad 7 ( t )) = 
Tr(ad^) with r\ = r y{t )~ 1 'j{to). Due to horizontality of 7, we have r] = Y^Li a iVh hence 
Tr (ad ,7) = Y%Li a tTr (ad p J = 0. Then the modular function is piecewise constant along 7. Re- 
calling that it is continuous, we have that it is constant. Varying along all g £ G and recalling 
that \l/(Id) = 1, we have \P = 1, hence G is unimodular. □ 



2.3.1 The hypoelliptic Laplacian on a non-unimodular Lie group 

In this section we present a non-unimodular Lie group endowed with a left-invariant sub- 
Riemannian structure. We then compute the explicit expression of the intrinsic hypoelliptic 
Laplacian: from Proposition 18 we already know that it is the sum of squares plus a first order 
term. 
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Consider the Lie group 



.4 H 



17 \ 

a b 
1 c 



iv 001 / 



a > 0, b, c e M > 



It is the direct sum of the group A + (M) of affine transformations on the real line x 
with a > and the additive group (R, +). Indeed, observe that 



ax + b 



( \ 

aOb 



1 c 



x 



( \ 



ax 



c + d 



v ooi nv i 1 / 



The group is non-unimodular, indeed a direct computation gives fiL 
-da db dc. 

a 

Its Lie algebra a(R) © M is generated by 



-^dadbdc and [ir 



Pi 



( \ 
1 





j 



P2 



( 

1 
1 



Oy 



/ \ 

1 




^0 Oy 



for which the following commutation rules hold: [pi,P2] — k [p2, k] = [k,pi] = —k. 
We define a trivializable sub-Riemannian structure on A + (IR) ©fas presented in Section 2.1.1: 
consider the two left-invariant vector fields Xi(g) = gpi with i — 1, 2 and define 
A(g) = span{Xi(#),X 2 (#)} g g (Xi(g),Xj(g)) = <J y . 



Using (9), one gets the following expression for the hypoelliptic Laplacian: 
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A sr = L 2 Xi (f) + L\ 2 <p + L Xl 4>. 

2.4 The intrinsic hypoelliptic Laplacian in terms of the formal adjoints of the vector fields 

In the literature another common definition of hypoelliptic Laplacian can be found (see for 
instance [30]): 

m 

a* = -£^*> (ii) 

i=l 

where {Xi, . . . , X m } is a set of vector fields satysfying the Hormander condition and the formal 
adjoint L* x . is computed with respect to a given volume form. This expression clearly simplifies 
to the sum of squares when the vector fields are formally skew-adjoint, i.e. L* x . = —L Xi - 
In this section we show that our definition of intrinsic hypoelliptic Laplacian coincides locally 
with (11), when {Xi, . . . , X m } is an orthonormal frame for the sub-Riemannian manifold and 
the formal adjoint of the vector fields are computed with respect to the sub-Riemannian volume 
form. 

We then show that left-invariant vector fields on a Lie group G are formally skew-adjoint with 
respect to the right-Haar measure, providing an alternative proof of the fact that for unimod- 
ular Lie Groups the intrinsic hypoelliptic Laplacian is the sum of squares. 

Proposition 19 Locally, the intrinsic hypoelliptic Laplacian A sr can be written as — Y%Li L* x .L Xi , 
where {Xi, . . . ,X m } is a local orthonormal frame, and L* Xj is the formal adjoint of the Lie 
derivative L Xi of the vector field X i; i.e. 

(01, = (02,^00, for every 1; 2 G C C °°(M, M), z = l,...,m, (12) 

and the scalar product is the one of L 2 (M, M) with respect to the invariant volume form, i.e. 

(01,02) := A/ 01 02 fJ-sr- 

Proof. Given a volume form \x on M, a definition of divergence of a smooth vector field X 
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(equivalent to L x /i = div(X)/i) is 

/ div(X)0/i = - / LxM for every E C C °°(M,R); 

JM JM 

see for instance [44]. We are going to prove that 

m 

A sr = - £ L^L Xi 0, for every E C C °°(M, R) : (13) 

i=l 

indeed, multiplying the left-hand side of (13) by if) E C^°(M) and integrating with respect to 
ji sr we have, 



/ (A sr 0) if; n„ = / (div sr (grad sr 0)) if) fi sr = / div sr YVLx^)^ if) /i 

JM JM JM \~[ J 

„ n 

Jm i=1 



For the right hand side we get the same expression. Since if) is arbitrary, the conclusion follows. 
Then, by density, one concludes that A sr (f> = - Y%Li L Xi L x .<fi, for every E C 2 (M, R). □ 
Proposition 20 Let G be a Lie group and X a left-invariant vector field on G. Then Lx is 
formally skew-adjoint with respect to the right-Haar measure. 

Proof. Let E C^°(M, R) and X = gp (p E L, g E G). Since X is left-invariant and fiR is 
right-invariant, we have 



J G (Lx<P)(9o)M9o) = j G j t 
~ dt 



<p(g e tp ) n R {g Q ) = ^ 



t=Q 



t=0 



G 



0) 



G 



t=0 

Hence, for every 1; 2 £ C£°(M, R) we have 

= / Lx(<pi4>2) Hr= I L x (<j>i) 02 At/? + / 0i (Lxfa) = (02, L x 4>i) + (4>i, L x 4> 2 ) 

JG JG JG 

and the conclusion follows. □ 

For unimodular groups we can assume fi sr = \xl = fiR (cfr. Remark 16) and left-invariant vector 
fields are formally skew-adjoint with respect to fi sr . This argument provides an alternative proof 
of the fact that on unimodular Lie groups the hypoelliptic Laplacian is the sum of squares. 
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3 The Generalized Fourier Transform on unimodular Lie groups 

Let / G L 1 (M, M): its Fourier transform is defined by the formula 

/(A) = / f(x)t "\/.r. 
jr 

If / G L 1 (R, R) n L 2 (R, R) then / G L 2 (R, R) and one has 

/ \f(x)\ 2 dx = [ |/(A)| 2 ^, 

JR JR. Z7T 

called Parseval or Plancherel equation. By density of L^R) n L 2 (R,R) in L 2 (R,R), this 
equation expresses the fact that the Fourier transform is an isometry between L 2 (R, R) and 
itself. Moreover, the following inversion formula holds: 

/(*) = f /(A)e^, 

JR 111 

where the equality is intended in the L 2 sense. It has been known from more than 50 years 
that the Fourier transform generalizes to a wide class of locally compact groups (see for in- 
stance [14,19,25,26,32,45]). Next we briefly present this generalization for groups satisfying the 
following hypothesis: 

(Ho) G is a unimodular Lie group of Type I. 
For the definition of groups of Type I see [18]. For our purposes it is sufficient to recall that all 
groups treated in this paper (i.e. H 2 , SU(2), 50(3), SL{2) and SE(2)) are of Type I. Actually, 
both the real connected semisimple and the real connected nilpotent Lie groups are of Type I 
[17,24] and even though not all solvable groups are of Type I, this is the case for the group 
of the rototranslations of the plane SE{2) [43]. In the following, the LP spaces L P (G,C) are 
intended with respect to the Haar measure \i := //£ = 

Let G be a Lie group satisfying (Ho) and G be the dual 6 of the group G, i.e. the set of all 
equivalence classes of unitary irreducible representations of G. Let A G G: in the following we 
6 See footnote 3. 
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indicate by X a choice of an irreducible representation in the class A. By definition, X is 
a map that to an element of G associates a unitary operator acting on a complex separable 
Hilbert space Ti, x : 

G -> U(H X ) 

X x : 

9 ^ %\g)- 

The index A for 7i x indicates that in general the Hilbert space can vary with A. 
Definition 21 Let G be a Lie group satisfying (Ho), and f G L l (G,C). The generalized (or 
noncommutative) Fourier transform (GFT) of f is the map (indicated in the following as f or 
J-(f)) that to each element of G associates the linear operator on 7i x : 

/(A) := JF(/) := / f{g)X\g- 1 )d l i. (14) 

Notice that since / is integrable and X x unitary, then /(A) is a bounded operator. 

Remark 22 / can be seen as an operator from f@ Ti to itself. We also use the notation 

/ =!g />) 

In general G is not a group and its structure can be quite complicated. In the case in which G 
is abelian then G is a group; if G is nilpotent then G has the structure of W 1 for some n; if G is 
compact then it is a Tannaka category (moreover, in this case each 7i x is finite dimensional). 
Under the hypothesis (H ) one can define on G a positive measure dP(X) (called the Plancherel 
measure) such that for every / G L 1 (G, C) fl L 2 {G, C) one has 

/ \f(g)Mg) = / >r(/(A) o f(\y)dP(X). 

By density of L 1 (G, C) fl L 2 (G, C) in L 2 (G, C), this formula expresses the fact that the GFT is 

e 

an isometry between L 2 (G, C) and Jq HS , the set of Hilbert-Schmidt operators with respect 
to the Plancherel measure. Moreover, it is obvious that: 

Proposition 23 Let G be a Lie group satisfying (Ho) and f G L 1 (G, C) fl L 2 (G, C). We have, 
for each g G G 
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f(g)= lTrif(\)oX x (g))dP(\). (15) 

JG 

where the equality is intended in the L? sense. 

It is immediate to verify that, given two functions fi, f 2 G L 1 (G, C) and defining their convo- 
lution as 

(A*/ 2 )(<?) = / fi(h)f 2 (h- 1 g)dh, (16) 

JG 

then the GFT maps the convolution into non-commutative product: 

Wi*/ 2 )(A) = / 2 (A)A(A). (17) 

Another important property is that if <5id(fi0 is the Dirac function at the identity over G, then 
5 Id (\) = Id H x. (18) 

In the following, a key role is played by the differential of the representation £ A , that is the 
map 



dX x : X ^ dX x (X) := 4 
v ; dt 



X\e tp ), (19) 

t=o 

where X = gp, (p G L, g G G) is a left-invariant vector field over G. By Stone theorem (see for 
instance [45, p. 6]) dX x (X) is a (possibly unbounded) skew-adjoint operator on Ti x . We have 
the following: 

Proposition 24 Let G be a Lie group satisfying (Ho) and X be a left-invariant vector field 
over G. The GFT of X , i.e. X = TLxT~ x splits into the Hilbert sum of operators X x , each 

one of which acts on the set HS A of Hilbert- Schmidt operators overH, x : 

e 

X = [ X x . 

JG 

Moreover, 

X X E = dX x (X) o S, for every E G HS A , (20) 
i.e. the GFT of a left- invariant vector field acts as a left-translation over HS A . 
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Proof. Consider the GFT of the operator R e t P of right-translation of a function by e tp , p G L, 
i.e. 

(iW) (9o) = f(9oe tp ), 

and compute its GFT: 

F(R etP f)(\)=F(f(g e t n) (A)= / /foe*)*^" 1 )^) = I /(^V^V^M*/^) 

J G j G 

= (#(e*)) /(A), 

where in the last equality we use the right-invariance of the Haar measure. Hence the GFT acts 
as a left-translation on HS A and it disintegrates the right-regular representations. It follows 
that 

© 

R e t P = J-R e t P J- 1 = / 3C x (e tp ). 

Jg 

Passing to the infinitesimal generators, with X = gp, the conclusion follows. □ 

Remark 25 From the fact that the GFT of a left-invariant vector field acts as a left-translation, 

it follows that X x can be interpreted as an operator over 7i x . 

3. 1 Computation of the kernel of the hypoelliptic heat equation 

In this section we provide a general method to compute the kernel of the hypoelliptic heat equa- 
tion on a left-invariant sub-Riemannian manifold (G, A, g) such that G satisfies the assumption 
(H ). 

We begin by recalling some existence results (for the semigroup of evolution and for the corre- 
sponding kernel) in the case of the sum of squares. We recall that for all the examples treated 
in this paper, the invariant hypoelliptic Laplacian is the sum of squares. 
Let G be a unimodular Lie group and (G, A,g) a left-invariant sub-Riemannian manifold 
generated by the orthonormal basis {pi, . . . ,p m }, and consider the hypoelliptic heat equation 
d t <f>{t,g) = A sr( f>{t,g). (21) 

Since G is unimodular, then A sr = L 2 Xi + . . . + L'j Cm , where Lx t is the Lie derivative w.r.t. 
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the vector field Xi := gpi (i — 1, . . . , m). Following Varopoulos [47, pp. 20-21, 106], since A sr 
is a sum of squares, then it is a symmetric operator that we identify with its Friedrichs (self- 
adjoint) extension, that is the infinitesimal generator of a (Markov) semigroup e tAsr . Thanks 
to the left-invariance of Xi (with i = 1, . . . , m), e tAsr admits a a right-convolution kernel Pt(-), 
i.e. 



„tA s 



M9) = 0o * Pt(g) = / Uh)pt{h- l g)^h) (22) 



G 



is the solution for t > to (21) with initial condition 0(0, g) = 0o(fiO £ L 1 (G, M) with respect 
to the Haar measure. 

Since the operator — A sr is hypoelliptic, then the kernel is a C°° function of (t, g) e M + x G. 
Notice that p t (g) = e tAsr Si d (g). 

The main results of the paper are based on the following key fact. 

Theorem 26 Let G be a Lie group satisfying (Ho) and (G, A, g) a left-invariant sub-Riemannian 
manifold generated by the orthonormal basis {pi, . . . ,p m }. Let A sr = L 2 Xi + . . . + L\ m be the 
intrinsic hypoelliptic Laplacian where L Xi is the Lie derivative w.r.t. the vector field Xi := gp^. 
Let {^ A } A ^ ^ e ^ e se t of all n on- equivalent classes of irreducible representations of the group 
G, each acting on an Hilbert space 7i x , and dP(X) be the Plancherel measure on the dual space 
G. We have the following: 

i) the GFT of A sr splits into the Hilbert sum of operators A^ r; each one of which leaves 7i x 
invariant: 

A sr = FA^F- 1 = f AidP(X), where A x sr = £ (x, A ) 2 . (23) 
Jg i=i 

ii) The operator A x r is self-adjoint and it is the infinitesimal generator of a contraction 
semi-group e tA ° r over HS A , i.e. e' A ^ r S A is the solution for t > to the operator equation 
d t E x (t) = A x r E x (t) in HS A , with initial condition S A (0) = S A . 

iii) The hypoelliptic heat kernel is 
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p t (g)= f^Tr(e tA -X x (g))dP(X), t > 0. (24) 

Proof. Following Varopoulos as above, and using Proposition 24, i) follows. Item ii) follows 
from the split (23) and from the fact that GFT is an isometry between L 2 (G,C) (the set 
of square integrable functions from G to C with respect to the Haar measure) and the set 

® A 

Jq HS of Hilbert-Schmidt operators with respect to the Plancherel measure. Item iii) is 
obtained applying the inverse GFT to e tA ^ r, E x and the convolution formula (17). The integral 
is convergent by the existence theorem for p t , see [47, p. 106]. □ 

Remark 27 As a consequence of Remark 25, it follows that A x r and e tA ° r can be considered 
as operators on 7i x . 

In the case when each A^ r has discrete spectrum, the following corollary gives an explicit 
formula for the hypoelliptic heat kernel in terms of its eigenvalues and eigenvectors. 
Corollary 28 Under the hypotheses of Theorem 26, if in addition we have that for every X, 
K x r (considered as an operator over7i x ) has discrete spectrum and |^n} ^ s a complete set of 
eigenf unctions of norm one with the corresponding set of eigenvalues jo^ j, then 

Ptis) = J 6 (E^n\2^)^>) dP(X) (25) 
where (., .) is the scalar product in 7i x . 

Proof. Recall that Tr (AB) = Tt(BA) and that Tr (A) = J2iei( e i, Aej) for any complete set 
{ e i)i & i °f orthonormal vectors. Hence Tr (^e t ^ r X x (g)^j = Y^ni^Pn' •£ A (flO e *^'"^n)- Observe that 
dtipn = ^sr^Pn = ^n^n) hence e tAsr i[) x = e a ™Vn5 from which the result follows. □ 



The following corollary gives a useful formula for the hypoelliptic heat kernel in the case in 
which for all A G G each operator e* A » r admits a convolution kernel Q x (., .). Here by ip x , we 
intend an element of 7i x . 

Corollary 29 Under the hypotheses of Theorem 26, if for all X £ G we have 7i x = L 2 (X X , d9 x ) 
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for some measure space (X , d6 ) and 



(9) = I i> x (9)Q x (9,9)d9 



then 



Pt(g) 



G JX X 



X\g)Q 



t V V 5 



d9dP(X), 



where in the last formula X x (g) acts on Q x (9,9) as a function of 9. 
Proof. From (24), we have 



p t (g) = Jjr {e^X x {g)) dP{\) = Jjr (X A (^)e^j dP{\). 



We have to compute the trace of the operator 



6 = X x (g)e m -- : ^ x (6) i-> X x (g)e m ^^ x (9) = X x (g) J ^ X (9)Q X (9, 9) d6 



x x 



K(9,9)^ x (9)d9 



(26) 



where K(9, 9) = X x (g)Q x (9, 9) is a function of 9, 9 and X x (g) acts on Q x (9, 9) as a function of 
9. The trace of G is f x K(9, 9)d9 and the conclusion follows. □ 



4 Explicit expressions on 3D unimodular Lie groups 

4-1 The hypoelliptic heat equation on H 2 

In this section we apply the method presented above to solve the hypoelliptic heat equation 
(21) on the Heisenberg group. This kernel, via the GFT, was first obtained by Hulanicki (see 
[29]). We present it as an application of Corollary 29, since in this case an expression for the 
kernel of the GFT of this equation is known. 
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We write the Heisenberg group as the 3D group of matrices 



1 x z + \xy 



Ho = 



1 y 



iv 00 1 



x,y,ze 



endowed with the standard matrix product. H 2 is indeed M 3 , 

/ \ 

1 x z + \xy 



(x,y,z) 



1 y 



v 00 1 / 



endowed with the group law 

(%i,yi, zi) ■ (x 2 , y 2 , z 2 ) = fxi + x 2 , yi + y 2 , z x + z 2 + ^ (x x y 2 - x 2 y x 
A basis of its Lie algebra is {pi,£>2? k} where 



Pi 



( \ 




( \ 




( \ 


1 









1 





P2 = 


1 


k = 





v° ) 




v° ) 




v° °y 



(27) 



They satisfy the following commutation rules: [pi,p2] = k, \pi,k] = [p 2 ,k] = 0, hence H 2 is a 
2-step nilpotent group. We define a left-invariant sub-Riemannian structure on H 2 as presented 
in Section 2.1.1: consider the two left-invariant vector fields Xi(g) = gpi with i — 1,2 and define 
k(g) = span{Xi(#),X 2 ( 5 )} g g (Xi(g), Xj{g)) = <J y . 



Writing the group H 2 in coordinates (x, y, z) on M 3 , we have the following expression for the 
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Lie derivatives of X\ and X 2 : 

L Xl =d x - |<9 Z , Lx 2 =d y + |<9 2 . 

The Heisenberg group is unimodular, hence the hypoelliptic Laplacian A sr is the sum of 
squares: 

A sr <f)= (L 2 Xi + L 2 X2 ) <f). (28) 

Remark 30 It is interesting to notice that all left-invariant sub-Riemannian structures that 
one can define on the Heisenberg group are isometric. 

In the next proposition we present the structure of the dual group of H 2 . For details and proofs 
see for instance [33]. 

Proposition 31 The dual space of H 2 is G = |X A | A G m} ; whe 



xere 



n -> n 

x\g) ■ 

whose domain isTi = L 2 (R, C), endowed with the standard product < ipi, ip 2 >'■= Jr " l l , i{9)ip 2 (9) d9 
where d9 is the Lebesgue measure. 

The Plancherel measure on G is dP(X) = ^d\, where d\ is the Lebesgue measure on R. 
Remark 32 Notice that in this example the domain 7i of the representation X x does not 
depend on A. 

4-1.1 The kernel of the hypoelliptic heat equation 

Consider the representation 3L X of H 2 and let X x be the corresponding representations of the 
differential operators L Xi with i = 1,2. Recall that X x are operators on Ti. From formula (19) 
we have 

[X^\{9) = ~m, [XfrW = -Mm, ^nce [A x sr i;}(9) = - A 2 ^ 2 ) 
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The GFT of the hypoelliptic heat equation is thus 



A 2 

d6 2 



The kernel of this equation is known (see for instance [4]) and it is called the Mehler kernel 
(its computation is very similar to the computation of the kernel for the harmonic oscillator 
in quantum mechanics): 

~ / 1 Acosh(2At) /q2 , ^ , \ee \ 



Q*( 6 >^ : \/27rsinh(2At) 6XP V 2 sinh(2At) v ^ ' " ' ' sinh(2At); ' 



+ e 2 ) + - 



Using Corollary 29 and after straightforward computations, one gets the kernel of the hypoel- 
liptic heat equation on the Heisenberg group: 



Pt{x,y,z) 



It 



exp 



r{x 2 + y 2 ) \ (0 zr 

— — - cos 2 — dr. 

2ttanh(2r) J K t 1 



(29) 



(2vrt) 2 Jr sinh(2r) 

This formula differs from the one by Gaveau [21] for some numerical factors since he studies 
the equation 

1 



dt<f> = \ [ip x + 2yd z ) 2 + (d y - 2xd z ) 2 
The Gaveau formula is recovered from (29) with t — > t/2 and z — > z/4. Moreover, a multiplica- 
tive factor | should be added, because from the change of variables one gets that the Haar 
measure is \dxdydz instead of dxdydz as used by Gaveau. 

4-2 The hypoelliptic heat equation on 577(2) 



In this section we solve the hypoelliptic heat equation (21) on the Lie group 



SU(2) 



{ ( \ 

a j3 



a,0e C, \a\ 2 + \(3\ 2 = 1 



A basis of the Lie algebra su(2) is {pi,P2, k) where 7 



7 See [43, pp. 67]: Pl = X u p 2 = X 2 , k = X 3 . 
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Pi 



( \ 

i 



v ! °/ 



P2 



( \ 

-1 



V 1 °J 



( \ 

i 



V°-'7 



(30) 



We define a trivializable sub-Riemannian structure on SU(2) as presented in Section 2.1.1: 
consider the two left-invariant vector fields Xi(g) = gpi with i = 1, 2 and define 

k(g) = span{X 1 (g),X 2 (g)} S g (Xi(g), Xj{g)) = <%. 



The group SU(2) is unimodular, hence the hypoelliptic Laplacian A sr has the following ex- 
pression: 



(31) 



In the next proposition we present the structure of the dual group of SU(2). For details and 
proofs see for instance [43]. 

Proposition 33 The dual space of 377(2) is G = {X n | n e N}. 

The domain 7i n of X n is the space of homogeneous polynomials of degree n in two variables 
(zi,Z2) with complex coefficients H n := {j2k=o a kZ k 
product 



k z 2 1 ~ k | a k e Cj, endowed with the scalar 



a k z\zT\ E Kz\zl~ k ) := E fc! (n - A;)! a k b k . 

k=0 k=0 k=0 



The representation 3£ n is 



x n (g) ■■ 



k jn—k 



k„.,n—k 



k=0 



E a k w k w 2 

k=0 



with (w 1 ,w 2 ) = (z 1 ,z 2 )g = (azi - (3z 2 ,(3z 1 + az 2 ). 

The Plancherel measure on G is dP{n) = (n + l)d^{n), where is the counting measure. 
Notice that an orthonormal basis of 7i n is {ibtVl n with ib2 := , 1 2 
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4-2.1 The kernel of the hypoelliptic heat equation 

Consider the representations X™ of the differential operators Lx t with i = 1,2: they are 
operators on 7i n , whose action on the basis {ipkYk=o °^ * s ( usrn g formula (19)) 

X^l = \ {Wk-i + (n- X^l = \ {k^U ~ (n - fc)Cfi} 



hence AJ.^ = [k 2 — kn — ~J ?p%. Thus, the basis {ipk}k=o is a complete set of eigenfunctions 
of norm one for the operator A™,. We are now able to compute the kernel of the hypoelliptic 
heat equation using formula (25). 

Proposition 34 The kernel of the hypoelliptic heat equation on (SU(2), A,g) is 



Pt 



where 



(32) 



n=0 



fc=0 



min{fc,n— k} 

A n > k (g):=(r k ,X n (g)r k )= E 

z=o 



k-l , 

V / 



n — k 



\ J 



S*-'a n -*-' lal 2 -l 



with g 



( \ 

a (3 



Proof. The formula p t {g) = T,^=o( n + 1) Efc=o e(fc ~ kn ~^ t {i } ki ^"(flO^jfc) is given by applying 
formula (25) in the SU(2) case. 

We now prove the explicit expression for (ip^,X n (g)ip^): a direct computation gives 



Y^min{s,n— fc} 
^--'Z=max{0,s— k} 



( \ 

k 



V s "'/ 



/ \ 

n — k 



\ 



a s-l^_^k-s+l^l^n-k-l 



\ J 



y/k\ {n-k)\ 
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Observe that ip^ is an orthonormal frame for the inner product: hence 



min{fc,n— k} 
1=0 



k-l , 
V / 



n — k 



V / 



a k -\-p) l (3 l a n - k ~ l ). 



The result easily follows. □ 

Remark 35 Notice that, as the sub-Riemannian distance (computed in [10]), pt(g) does not de- 
pend on (3. This is due to the cylindrical symmetry of the distribution around e k = |e cfc | c e M j. 



4-3 The hypoelliptic heat equation on 50(3) 

Let g be an element of 50(3) = [A e Mat(R,3) | AA T = Id, det(A) 
Lie algebra so(3) is {pi,P2, k} where 8 



l|. A basis of the 



Pi 



( \ 


0-1 



v 01 °j 



Pi 



( \ 
1 





V" 100 / 



0-10 
1 



Oj 



(33) 



We define a trivializable sub-Riemannian structure on 50(3) as presented in Section 2.1.1: 
consider the two left-invariant vector fields Xi(g) = gpi with % = 1, 2 and define 
k(g) = span{X 1 (^),X 2 ( 5 )} g g (Xi(g), Xj{g)) = <J y . 

The group 50(3) is unimodular, hence the hypoelliptic Laplacian A sr has the following ex- 
pression: 



L 2 Xl + L 2 x U. (34) 



We present now the structure of the dual group of 50(3). For details and proofs see [43]. 



See [43, pp. 88]: pi = Z\, p 2 = Z 2 , k = Z 3 . 
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First consider the domain Ti. r , that is the space of complex- valued polynomials of r-th degree 
in three real variables x, y, z that are homogeneous and harmonic 
7f' = {f(x, y, z) | deg (/) = r, f homogeneous , Af = 0} . 

Notice that an homogeneous polynomial / G 7i r is uniquely determined by its value on S 2 = 

{(x,y,z) | x 2 + y 2 + z 2 = 1}, as f(px,py,pz) = p r f(x,y,z). 

Define f(a,/3) := /(sin(a) cos(/5), sin(a) sin(/3), cos(a)). Then endow 7i r with the scalar prod- 
uct 

< h(x,y,z),f 2 (x,y,z) >:= — f f 1 (a,p)f 2 (a,f3) sin a da d/3. 

47T JS 2 

In the following proposition we present the structure of the dual group. 
Proposition 36 The dual space of SO (3) is G = {X r \ r G N}. 

Given g G 5*0(3), the unitary representation 3L r (g) is 

x r (g) ■ 

f(x,y,z) i * f(xi,Vh *i) 

with (si,2/i,^i) = (x,y,z)g. 

The Plancherel measure on 50(3) zs dP(r) = (2r+l)dp,$(r), where dp$ is the counting measure. 
An orthonormal basis for H r is given by {(p r s } r s= _ r with <$.(a, (3) := e lsl3 P~ s (cos(a)), where 
P*(x) are the Legendre polynomials 9 . 

4-3.1 The kernel of the hypoelliptic heat equation 

Consider the representations X\ of the differential operators with % — 1,2: using formula 

(19) we find the following expressions in spherical coordinates 10 

dip dip " dip dip 

X[ip = sin(/5)— + cot (a) cos(/3)— , X r 2 ip = - cos(/5)— + cot(a) sin(/3) — 

hence 

9 Recall that P r s (x) is defined by P r s (x) := ^7^- d + j£r+l^ ■ 

10 i.e. x = psin(a) cos(/3), y = psin(a) sin(/3), z = pcos(a). 
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A W = g +00t V,g + (35) 

and its action on the basis {4> r s } r s= _ r of W is 

A^=( S 2 -r(r + l))^. (36) 

Hence the basis {0s}^ = _ r is a complete set of eigenfunctions of norm one for the operator A" r . 
We compute the kernel of the hypoelliptic heat equation, using (25). 
Proposition 37 The kernel of the hypoelliptic heat equation on (50(3), A, g) is 

oo r 

p t (g) = £(2r + 1) £ e(* 2 -^+D)' < fl, X r (g)<f>: > . (37) 

r=0 s=—r 

4-3.2 The heat kernel on 50(3) via the heat kernel on SU(2) 

In this section we verify that the heat kernel on 50(3) given in (37) can be easily retrieved 
from the one on SU(2) given in (32). In the following, all the objects relative to 50(3) are 
underlined, e.g. g G 50(3),p, G so(3), the representations XT. acting on JV_ with basis (j) r s . 
Consider the isomorphism of Lie algebras ad : su(2) — > so(3) defined by adpi = pi, ad^2 = 
P2j ad/c = fc: it gives the matrix expression of the adjoint map on su{2) with respect to the 
basis {pi,p2,k}. There is a corresponding endomorphism of groups Ad : SU{2) — > 50(3) 
given by Ad(exp-u) = exp(ad(t>)). It is a covering map of 50(3) by SU(2), such that for each 
matrix g G 50(3) the preimage is given by two opposite matrices g, —g G SU{2). 
Proposition 38 The following relation holds between the kernel pt on 50(3) given in (37) 
and the kernel p t on SU(2) given in (32): 
\/geS0(3),geAd~\g) 

, , Pt(g) +p t (-g) 
Pt{g) = 2 • 

Proof. Observe the following key facts (see e.g. [43, II. §7]): 

• on SU{2): ^(-g)<f> = (-l) n X n (g)fr 

• the representation 37 of 50(3) and the representation X 2r of SU(2) are unitarily related, 
i.e. the following relation holds: V g G SU{2), (jf s G JV_ 
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(38) 



where the map 



n 



'li- 



ar 
'r+s- 



is a unitary isomorphism. 
Then we have explicitly 

Pt(g)+Pt(-9) _ S£ (n +!)(! + (-!)") E n k=0 e^- kn -^A n ' k (9) _ 
2 2 

oo r 

= £(2r + 1) E e (s2 - r(r+1))t < 'f&.t&Wtfc. > 

r=0 s=—r 

after the substitution r = 2n, s = A; — r. Using Equation (38), we have < 
(f) r s , X[_(Adg)<f) r s >, from which the result directly follows. □ 



i2r Tlr( „\a3.t 



x 2r (gW; +s >=< 



4-4 The hypoelliptic heat equation on SL(2) 



In this section we solve the hypoelliptic heat equation (21) on the Lie group 



SL(2) = {ge Mat(R, 2) | det(g) = 1} 



A basis of the Lie algebra si (2) is 

/ 

1 



Pl= 2 



0-1 



P2= 2 



( \ 
1 



v 10 / 



2 



-1 



v 1 V 



We define a trivializable sub-Riemannian structure on SL{2) as presented in Section 2.1.1: 
consider the two left-invariant vector fields X^g) = gpi with i — 1,2 and define 
A(#) = spm{X 1 (g),X 2 (g)} g g (Xi(g),Xj(g)) = 6 tj . 



The group SL(2) is unimodular, hence the hypoelliptic Laplacian A sr has the following ex- 
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pression: 



A sr 0= (L 2 Xl +L 2 Xa ) 



(39) 



It is well known that SL(2) and 577(1, 1] 
groups via the isomorphism 



if \ 

a (3 



a 



1 > are isomorph Lie 



SL(2) 



SU {1,1) 



IT : 



with C 



g ^ Sf = CgC" 1 



V2 



( \ 
1 -i 



v 1 1 / 



This isomorphism also induce an isomorphism of Lie algebras dll : sl(2) — > su(l, 1) defined by 
effl(pi) = pi, dn(p 2 ) = j/ 2 , rfn(fc) = k' with 



1 



v 10 / 



P2= 2 



( \ 

-i 



k' = - 
2 



( \ 

-i 



This isomorphism induces naturally the definitions of left-invariant sub-Riemannian structure 
and of the hypoelliptic Laplacian on SU(1, 1). 

We present here the structure of the dual of the group SU(1, 1), observing that the isomorphism 

of groups induces an isomorphism of representations. For details and proofs, see [43] . 

The dual space G of SU(1, 1) contains two continuous and two discrete parts: G = Gq II Gu 

with G c = {&> s | j E {0, 1} , s = \ + iv,v e M + } and G D = [x n \ n e §Z, |n| > l}. 

We define the domain Tic of the continuous representation j£- ?,s : it is the Hilbert space of 

L 2 complex- valued functions on the unitary circle S 1 = {i 6 C |x| = l} with respect to 

the normalized Lebesgue measure endowed with the standard scalar product (/, g) : = 



Is 1 f( x )d( x ) |f- A n orthonormal basis is given by the set {^ m } me z with <jf s 



x) 



X 
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Proposition 39 The continuous part of the dual space of SU(1, 1) is 



G 



o 



Given & e SU(1, 1), the unitary representation X 



3,s 



is 



n 



c 



n 



c 



\(3x + a\ 



-2s 



with C S' 



I \ 

a (3 



fix + a 
\/3x + a\ 



2j 



'ax + /5 N 
fix + a 



The Plancherel measure on Gc is dP(j, \ + iv) 



^Tanh(7Tf ) dv j ' = 
^:Cotanh(7r?;) dv j = 



where dv 



is the Lebesgue measure on R. 

Remark 40 Notice that in this example the domain of the representation Tic does not depend 
on j, s. 

Now we turn our attention to the description of principal discrete representations 11 . 

We first define the domain 7i n of these representations X n : consider the space S£ n of L 2 

complex- valued functions on the unitary disc I> = {i G C |a;| < 1} with respect to the mea- 



sure dfj,*(z) 



2 n -1 , 



|2\2n-2 



dz where dz is the Lebesgue measure. ££ n is an Hilbert space 



if endowed with the scalar product (/, g) := f D f(z)g(z) d/j,*(z). Then define the space H n with 
n > as the Hilbert space of holomorphic functions of J?f n , while 7i n with n < is the Hilbert 
space of antiholomorphic functions of J2?_ n - An orthonormal basis for 7i n with n > is given 
by {OmgN with ip^iz) = ( r(2n)r(m+i) ) 2 zm w h er e r is the Gamma function. An orthonormal 



11 There exist also the so-called complementary discrete representations, whose Plancherel measure is 
vanishing. Hence they do not contribute to the GFT of a function defined on SU(1, 1). For details, 
see for instance [43]. 
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basis for H n with n < is given by {^m}meN w ith = ip-n( z )- 

Proposition 41 The discrete part of the dual space ofSU(l, 1) is Gn = |X n | n G |Z, \n\ > l| 
Given £f G SU(1, 1), £/ie unitary representation £ n (£f) is 



-1 



a /3 



The Plancherel measure on Go is dP{n) = 1 d/j,$(n), where d/xj i/ie counting measure. 
4-4-1 The kernel of the hypoelliptic heat equation 

In this section we compute the representation of differential operators Lx t with % = 1,2 and 
give the explicit expression of the kernel of the hypoelliptic heat equation. 
We first study the continuous representations Xj' s , for both the families j = 0, |. Their actions 
on the basis {^ m } mgZ of He is 

■tris, s — m — j. s + m + j 

X l Wm = ~ Wrn-1 H Z Wm+1, 



yj, s i .s-m-j .s + m + j 



Hence 



A j,s tb 



-\[m + j, 



Moreover, the set {?pm} m( zz is a complete set of eigenfunctions of norm one for the operator 
A j ' s . 

sr 

Remark 42 Notice that the operators Xf' s are only defined on the space of C°° vectors, i.e. 
the vectors v G Tic such that the map g — > [X^ ,s (g)} v is a C°° mapping. This restriction is not 
crucial for the following treatment. 

We now turn our attention to the discrete representations in both cases n > (holomorphic 
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functions) and n < (antiholomorphic functions). Consider the discrete representation X n of 
SU(1, 1) and let X? be the representations of the differential operators Lx t with i = 1,2. Their 
actions on the basis {^} m6N of H n are 



\/(2|n| + m)(m + 1) </(2|n| + m — l)m 

A irra" o fro+1 o rm-1 



. yjj +m)(m + 1) V(2|n| +m- l)m 

z o rm+i — 1 r, rm-i 



Hence A™^™ = — (|n| + 2m|n| + m 2 ) thus the basis {ft} m6N is a complete set of eigen- 
functions of norm one for the operator A" r . 

We now compute the kernel of the hypoelliptic heat equation using Formula (25). 
Proposition 43 The kernel of the hypoelliptic heat equation on (SL(2), A,g) is 



Pt 



(<?) = r° ^TanhM £ e < m2+ ^) (rj; m , X°' s (^ m ) dfi(v) + 
+ f °°^CotanhM ]T e< m ^ mW ^\^ m ^ s {^ m ) d/i(v) + 



+ E E e"*('"i +2m i"' +m2 )(c^ n (^)C)- 



neiZ, |n|>l 



mGN 



w/iere^ = n^- 1 ) e Sf/(1,1). 



4-5 The hypoelliptic heat kernel on SE(2) 



(40) 



Consider the group of rototranslations of the plane 

it \ 

cos(a) — sin (a) x\ 



SE(2) 



sin « cost a) x<i 



o 







a e R/2%, Xi e 



J 



In the following we often denote an element of SE(2) as g = (a,Xi,X2) 
A basis of the Lie algebra of SE{2) is {po,Pi,P2}, with 
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Po 



/ \ 

0-10 
1 



v° °/ 



Pi 



/ \ 
1 





v / 



, P2 



/ \ 


1 



v / 



(41) 



We define a trivializable sub-Riemannian structure on SE(2) as presented in Section 2.1.1: 
consider the two left-invariant vector fields Xi(g) = gpi with i = 0, 1 and define 

k{g) = span {X (g), X^g)} g g (Xi(g), Xj(g)) = c%. 

The group SE(2) is unimodular, hence the hypoelliptic Laplacian A sr has the following ex- 
pression: 



A, 



L X + L Xx 



(42) 



Remark 44 As for the Heisenberg group, all left-invariant sub-Riemannian structures that one 
can define on SE(2) are isometric. 

In the following proposition we present the structure of the dual of SE(2). 
Proposition 45 The dual space of SE(2) is G = |X A | A G R+}. 

Given g = (a,xi,X2) G SE{2), the unitary representation £ A (g) is 



n 



n 



ip(8) ^ e iX(xcos( - e) - ysin( - e)) i:(9 + a) 



where the domain TL of the representation % x (g) is 71 = L 2 (S' 1 ,C) ; the Hilbert space of L 2 
functions on the circle S l C IR 2 with respect to the Lebesgue measure d9, endowed with the 



scalar product < ipi,ip2 >= Is 1 "^l (69^2(0) d6. 

The Plancherel measure on G is dP(X) = XdX where dX is the Lebesgue measure on R. 
Remark 46 Notice that in this example the domain of the representation 7i does not depend 
on A. 
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4-5.1 The kernel of the hypoelliptic heat equation 

Consider the representations X x of the differential operators Lx i with i = 1,2: they are ope- 
rators on 7i, whose action on ip £ 7i is (using formula (19)) 

(9) = *^p. [jtfr] (9) = i\ caB{0)m, 

hence 

\kr^\ W = ^|P-A 2 COS 2 (^). 



We have to find a complete set of eigenf unctions of norm one for A" r . An eigenfunction ip with 
eigenvalue E is a 27r-periodic function satisfying the Mathieu's equation 



dx 2 



+ (a-2gcos(2x))^ = (43) 



with a = — 4r — E and q = For details about Mathieu functions see for instance [1, ch. 20]. 
Remark 47 Notice that we consider only 27r-periodic solutions of (43) since TC = L 2 (S' 1 , C). 
There exists an ordered discrete set {a n (q)} n = of distinct real numbers (a n < a n +i) such that 
the equation + (a n — 2gcos(2x))/ = admits a unique even 27r-periodic solution of norm 
one. This function ce n (x,q) is called an even Mathieu function. 

Similarly, there exists an ordered discrete set {b n (q)} n =i of distinct real numbers (b n < b n+ \) 
such that the equation jjj + (b n — 2qcos(2x))f = admits a unique odd 27r-periodic solution 
of norm 1. This function se n (x, q) is called an odd Mathieu function. 

The set B x := |ce n (x, ^ | o Ll|se n (x, ^ | is a complete set of 27r-periodic eigenf unctions 
of norm one for the operator A" r . The eigenvalue for ce n (x, is a A := — 4p — a n (^-j . The 
eigenvalue for se n (x, ^ ) is 6 A := — ^ — fo r ^ 



' 4 / " ' 2 n V 4 

We can now compute the explicit expression of the hypoelliptic kernel on SE{2). 
Proposition 48 The kernel of the hypoelliptic heat equation on (SE(2), g) is 

hoo +oo 

~ A d\ ( ' 

'o 



/■ + 00 / + oo +oo \ 

Pt(^) = y A dA ( 5> a «* < ce n (0),X A (<7)ce n (0) > + ]T < se n (0), £ A Q7)se(fl) > (44) 



\n=0 n=l 
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The function (44) is real for allt > 0. 

Proof. The formula (44) is given by writing the formula (25) in the SE(2) case. 

We have to prove that pt{g) is real: we claim that < ce n , X x (g)ce n > is real. In fact, write the 

scalar product with g = (a,x,y): 

< ce n , X\g)ce n >= / e iX( - xcos ^ ys ^ce n (9)ce n (9 + a). 
Jo 

Its imaginary part is sin (A(x cos(6 l ) — y sin(#))) ce n (6)ce n (9 + a). Its integrand function 
assumes opposite values in 9 and 9 + n: indeed 

sin (X(x cos(6 ) + ir) — y sin(6' + 7r))) = sin (A(— x cos(6 ) ) + y sm(9))) = — sin (A(+x cos(6 ) ) — y sin(6 1 ))) , 

while ce n (# + tt) = (—l) n ce n (9) as a property of Mathieu functions. Thus, the integral over 
[0,27r] is null. With similar observations it is possible to prove that < se n (9), % x (g)se n (9) > is 
real. 

Thus, pt{g) is an integral of a sum of real functions, hence it is real. □ 

Acknowledgments. The authors are grateful to Giovanna Citti and to Fulvio Ricci for helpful 
discussions. 



15 



References 

[I] M. Abramowitz, I. A. Stegun (Editors), Handbook of Mathematical Functions with 
Formulas, Graphs and Mathematical Tables, National Bureau of Standards Applied Mathematics, 
Washington, 1964. 

[2] A. A. Agrachev, Yu.L. Sachkov, Control Theory from the Geometric Viewpoint, Encyclopedia 
of Mathematical Sciences, 87, Springer, 2004. 

[3] A. A. Agrachev, U. Boscain, M. Sigalotti, A Gauss-Bonnet-like Formula on Two- 
Dimensional Almost-Riemannian Manifolds, Discrete Contin. Dyn. Syst. 20, no. 4, 801-822, 2008. 

[4] R. Beals, Solutions fondamentales exactes. (French) [Exact kernels] Journees "Equations aux 
Derivees Partielles" (Saint- Jean-de-Monts, 1998), Exp. No. I, 9 pp., Univ. Nantes, Nantes, 1998. 

[5] R. Beals, B. Gaveau, P. Greiner, Hamilton-Jacobi theory and the heat kernel on Heisenberg 
groups. J. Math. Pures Appl. (9) 79, no. 7, pp. 633-689, 2000. 

[6] A. Bellaiche, The tangent space in sub-Riemannian geometry, in Sub-Riemannian geometry, 
edited by A. Bellaiche and J. -J. Risler, pp. 1-78, Progr. Math., 144, Birkhauser, Basel, 1996. 

[7] R. Brockett, Nonlinear control theory and differential geometry. Proceedings of the 
International Congress of Mathematicians, Vol. 1, 2, Warsaw, (1984), pp. 1357-1368. 

[8] R. Brockett, Control theory and singular Riemannian geometry, in New directions in applied 
mathematics (Cleveland, Ohio, 1980), pp. 11-27, Springer, New York-Berlin, 1982. 

[9] A. Bonfiglioli, E. Lanconelli, F. Uguzzoni, Stratified Lie groups and potential theory for 
their sub-Laplacians, Springer Monographs in Mathematics, Springer Berlin, 2007. 

[10] U. Boscain, F. Rossi, Invariant Carnot-Caratheodory metrics on S 3 , <SO(3), SL(2) and Lens 
Spaces, SIAM J. on Contr. and Optim., SIAM J. Control Optim., 47 (2008), no. 4, pp. 1851-1878. 

[II] U. Boscain, S. Polidoro, Gaussian estimates for hypoelliptic operators via optimal control, 
Atti Accad. Naz. Lincei CI. Sci. Fis. Mat. Natur. Rend. Lincei (9) Mat. Appl, 18 (2007), no. 4, 



16 



pp. 333-342. 

[12] A. Bressan, B. Piccoli, Introduction to the Mathematical Theory of Control, AIMS Series 
on Applied Mathematics n. 2, 2007. 

[13] L. Capogna, D. Danielli, S.D. Pauls, J.T. Tyson, An introduction to the Heisenb erg group 
and the sub-Riemannian isoperimetric problem. Progress in Mathematics, 259. Birkhuser Verlag, 
Basel, 2007. 

[14] G. S. Chirikjian, A. B. Kyatkin, Engineering applications of noncommutative harmonic 
analysis, CRC Press, Boca Raton, FL, 2001. 

[15] G. Citti, A. Sarti, A cortical based model of perceptual completion in the roto-translation 
space, AMS Acta, 2004. 

[16] J. Cygan, Heat kernels for class 2 nilpotent groups, Studia Math., 64 (1979), no. 3, pp. 227-238. 

[17] J. Dixmier, Sur les representations unitaires des groupes de Lie algebriques. Ann. Inst. Fourier, 
Grenoble, 7 (1957), pp. 315-328. 

[18] J. Dixmier, Sur les representations unitaires des groupes de Lie nilpotents. I. Amer. J. Math. 
81, (1959) pp. 160-170. 

[19] M. Duflo, Analyse harmonique sur les groupes algebriques complexes: formule de Plancherel 
(d'apres M. Andler) et conjecture de M. Vergne, Bourbaki seminar, Vol. 1982/83, pp. 279-291. 

[20] G.B. Folland, E. M. Stein, Estimates for the db complex and analysis on the Heisenberg 
group. Comm. Pure Appl. Math. 27 (1974), pp. 429-522. 

[21] B. Gaveau, Principe de moindre action, propagation de la chaleur et estimees sous elliptiques 
sur certains groupes nilpotents, Acta Math., 139 (1977), no. 1-2, pp. 95-153. 

[22] V. Gershkovich, A. Vershik, Nonholonomic manifolds and nilpotent analysis. J. Geom. Phys. 
5 (1988), pp. 407-452. 



17 



[23] M. Gromov, Carnot-Caratheodory spaces seen from within, in Sub-Riemannian geometry, edited 
by A. Bellaiche and J. -J. Risler, Progr. Math., 144, Birkhauser, Basel, 1996, pp. 79-323. 

[24] Harish-Chandra, Representations of a semisimple Lie group on a Banach space. I. Trans. Amer. 
Math. Soc. 75, (1953) pp. 185-243. 

[25] E. Hewitt, K.A. Ross, Abstract harmonic analysis. Vol. I. Structure of topological groups, 
integration theory, group representations. Second edition. Grundlehren der Mathematischen 
Wissenschaften, 115. Springer- Verlag, Berlin-New York, 1979. 

[26] E. Hewitt, K.A. Ross, Abstract harmonic analysis. Vol. II: Structure and analysis for compact 
groups. Analysis on locally compact Abelian groups. Die Grundlehren der mathematischen 
Wissenschaften, Band 152 Springer- Verlag, New York-Berlin, 1970. 

[27] H. Heyer, Dualitat lokalkompakter Gruppen. Lecture Notes in Mathematics, Vol. 150 Springer- 
Verlag, Berlin-New York 1970. 

[28] L. Hormander, Hypoelliptic Second Order Differential Equations, Acta Math., 119 (1967), 
pp. 147-171. 

[29] A. Hulanicki, The distribution of energy in the Brownian motion in the Gaussian field and 
analytic-hypoellipticity of certain subelliptic operators on the Heisenberg group, Studia Math. 56 
(1976), no. 2, pp. 165-173. 

[30] D. Jerison, A. Sanchez-Calle Subelliptic, second order differential operators, in Complex 
analysis, III, Lecture Notes in Math., 1277, Springer, Berlin, 1987, pp. 46-77. 

[31] D. Jerison, A. Sanchez-Calle, Estimates for the heat kernel for the sum of squares of vector 
fields, Indiana Univ. Math. Journ., 35 (1986), pp. 835-854. 

[32] A. A. Kirillov, Elements of the theory of representations. Grundlehren der Mathematischen 
Wissenschaften, Band 220. Springer- Verlag, Berlin-New York, 1976. 

[33] A. A. Kirillov, Lectures on the orbit method, American Mathematical Society 2004 Series 
Graduate studies in mathematics 64, 2004. 



48 



[34] A. Klingler, New derivation of the Heisenberg kernel Comm. Partial Differential Equations 22, 
1997, no. 11-12, pp. 2051-2060. 

[35] R. Leandre, Minoration en temps petit de la densit d'une diffusion dgnre. J. Funct. Anal. 74 
(1987), pp. 399-414. 

[36] J. Mitchell, On Carnot-Carathodory metrics. J. Differential Geom. 21 (1985), no. 1, pp. 35-45. 

[37] I. Moiseev, Yu.L. Sachkov, Maxwell strata in sub-Riemannian problem on the group of 
motions of a plane, arXiv:0807.4731vl. 

[38] R. Montgomery, A Tour of Subriemannian Geometries, their Geodesies and Applications, 
Mathematical Surveys and Monographs 91, AMS, Providence, RI, 2002. 

[39] R. Neel, D. Stroock, Analysis of the cut locus via the heat kernel, in Surv. Differ. Geom., 
IX, Int. Press, Somerville, MA, (2004), pp. 337-349. 

[40] P. Pansu Metriques de Carnot-Caratheodory et quasiisometries des espaces symetriques de rang 
un, Ann. of Math. (2) 129 (1989), no. 1, pp. 1-60. 

[41] J. Petitot, Vers une Neuro-geometrie. Fibrations corticales, structures de contact et contours 
subjectifs modaux, Numero special de Mathematiques, Informatique et Sciences Humaines, n. 
145, EHESS, Paris, 1999, pp. 5-101. 

[42] L. P. Rothschild, E. M. Stein, Hypoelliptic differential operators and nilpotent groups. Acta 
Math. 137 (1976), no. 3-4, pp. 247-320. 

[43] M. SUGIURA, Unitary representations and harmonic analysis: an introduction, 2nd Edition, 
North-Holland Mathematical Library 44, North-Holland Publishing Co., Amsterdam; Kodansha, 
Ltd., Tokyo, 1990. 

[44] M. E. Taylor, Partial differential equations. I. Basic theory. Applied Mathematical Sciences, 
115. Springer- Verlag, New York, 1996. 

[45] M. E. Taylor, Noncommutative harmonic analysis. Mathematical Surveys and Monographs, 
22. American Mathematical Society, Providence, RI, 1986. 



49 



[46] T. J. S. Taylor, Off diagonal asymptotics of hypoelliptic diffusion equations and singular 
Riemannian geometry. Pacific J. Math. 136 (1989), pp. 379-399. 

[47] N. VAROPOULOS, L. Saloff-Coste, T. Coulhon, Analysis and geometry on groups, 
Cambridge Tracts in Mathematics, 100, Cambridge University Press, Cambridge, 1992. 



50 



